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Abstract 

We re-consider the time dependent Schrodinger-Newton equation as a model for 
the self-gravitational interaction of a quantum system. We numerically locate 
the onset of gravitationally induced inhibitions of dispersion of Gaussian wave 
packets and find them to occur at mass values more than 6 orders of magnitude 
higher than reported by Salzman and Carlip [9, 2|, namely at about 10 10 u. This 
fits much better to simple analytical estimates but unfortunately also questions 
the experimental realisability of the proposed laboratory test of quantum gravity 
in the foreseeable future, not just because of large masses, but also because of 
the need to provide sufficiently long coherence times. 



1 Introduction 

In this paper we investigate the time dependent Schrodinger-Newton equation (hence- 
forth abbreviated SN-equation): 

ihdt*(t,x) = (-— A-Gm 2 [ ^i*'^|, 2 d 3 y) V(t,x). (1) 
V 2m J \\x-y\\ J 

It can be thought of as ordinary time-dependent Schrodinger equation, 

ihd t ^(t,x) = (-— A + m$(t,x)) (2) 
\ 2m J 

with a Newtonian gravitational potential <I> sourced by a mass distribution associated 
to \& itself, that is m|^| 2 , so that 

A$(t,£) = AnGm\^(t,x)\ 2 . (3) 
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The re-coupling of ^ via ([3]) into ([2]) results in the non-linear and non-local self- 
interaction seen in ([T]). Physical intuition suggests this self-interaction to cause a slow 
down of dispersion which should be the more pronounced the higher the value of m is 
chosen. 

In 10 it was surprisingly found on the basis of numerical computations that 
Gaussian wave packets of width a « 0.5 fim. start shrinking in width if the mass pa- 
rameter exceeds m m 1600 u. This is surprising in view of the fact that simple ana- 
lytical estimates suggest this type of behaviour to occur much later, namely at the said 
initial width for mass values between around 10 10 u. The authors of []9l |2| observed 
this discrepancy and encouraged a check of their results, which we performed. The 
result is, that our numerical studies now fully confirm the high mass values so that no 
discrepancy between numerical and simple analytical estimates seem to occur. 

The special attention to typical width dimension of a « 0.5 fj,m stems from actual 
molecular interferometry experiments in which the wave nature of complex molecules 
has been demonstrated; see, e. g., HI for an overview. The molecules used in [4] were 
tetraphenylporphyrin C^H^qN^ with mass 614 u and the more complex and more 
massive fluorofullerene CqqF^ composed of 108 atoms and of total mass 1632 u. The 
grating period and grating thickness were about 991 nm and 500 nm respectively. 

For further discussion, we now collect some background material concerning the 
SN-equation ([T]). 

2 Basic features of the Schrodinger-Newton Equation 

Equation {!]) can be derived from the following action (we write & instead of and 
^* for the complex conjugate of ^): 

$[¥,¥*] = J dijy j d 3 x(V(i,x)^(t,f) -V(t,2)4f*(t,x 

/ d 3 x(V^(t,x)) ■ (V^*(t,x)) 

2m J v ' 

o™ 2 rrM-j>. m,z)\ 2 m -•" 



Note that VI/ and \&* are to be varied independently. A local form can be obtained if the 
Newtonian gravitational field is introduced: 

**,*] = j dt J d 3 x!^(^*(t,x)if(t,x) - V(t,x)V*(t,x) 

-^(^(t,x)).(V^(t,x)) 

1 ;{V<5>(t,x)) • (V#(i,x))) 



8ttG 

m\^(t,x)\ 2 <S>(t,x) } . (5) 



Variation of with respect to gives Inserting the solution = A-KGmA 1 1 \? | 2 
into ([5) then leads to Q. 
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We usually restrict attention to solutions which are square integrable over space, 
so that solutions of the time dependent SN-equations are paths in the Hilbert space 
L 2 (M 3 ,d 3 x). We shall see below that time evolution preserves the norm (this follows 
from the symmetry under phase transformations). 

Symmetries of ([T]) are: 

1 . Constant phase shifts 

* i->- , H?'(t, x) := exp(ia)V(t, x) (6) 

with constant a G K. This symmetry of the equation of motion is also a sym- 
metry of the action ([4]). Hence there is a conserved Noether current which turns 
out to be just the same functional of ^ as in ordinary Schrodinger theory. In 
particular, the space integral over |^| 2 is time independent. 

2. Proper orthochroneous Galilei Transformations 

(t, x) -)• T g (t, x ) := (t + b , R • x + vt + a) , (7) 

where R G SO(3) denotes spatial rotations, v G M 3 boost transformations, 
and a G M 3 and b G M. space and time translations respectively, and where g 
collectively stands for (R, v, a, b). These generate the inhomogeneous (proper 
orthochronous) Galilei group which acts via proper ray -representations on wave 
functions as follows: 

* -> T g ^ := exp(i/3 9 )(^ o Tg -x) (8) 

Note that except for the phase exp(i/3 g ), this is just the transformation of a scalar 
function on spacetime under the diffeomorphism r g . The (7-dependent phase 
function (3 : M 4 — > E is determined up to spacetime independent additions. A 
convenient choice is (compare Q) 



f3 g (t,x) = r ^ v-{x-a)-\v 2 {t-b) 



(9) 



It is immediate that the self-interaction term in ([T]) and Q does not affect Galilei 
invariance of the equation of motion. Moreover, ([8]) is also a symmetry of the 
action so that 10 conserved quantities result. The expressions for total momen- 
tum and angular momentum are again the same functionals of \£ as in ordinary 
Schrodinger theory, but that for total energy will clearly include the non-linear 
interaction. 

3. Whereas the proper orthochroneous Galilei transformations exclude space and 
time inversions, the latter are also symmetries of the SN-equation. To define 
them, let ir and r be the following diffeomorphisms of M 4 : 

7r : M 4 — t- M 4 , (t, x) ^ (t, -x) , (10a) 
r:M 4 ^]R 4 , (t, x) ^ (-t, x) . (10b) 
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The unitary transformations of space inversion and the anti-unitary transforma- 
tion of time inversion (or better "reversal of motion") on wavefunctions is then 
defined by 

Ptf := ¥ o tt -1 , (11a) 
Tf ^Cofor 1 , (lib) 

where C denotes complex-conjugation. It is easy to check that if ^ solves the 
SN-equation then Pty and also solve it. 

There is a scaling equivariance between SN-equations of different mass param- 
eters m. The precise statement is this: Let IR + be the multiplicative group of 
positive real numbers. It acts by scaling transformations on spacetime: 



<T '. K+ X K — > K , a 



(// , (t, x)) = a^t, x) = (/i a t , n h x) , (12) 



where a, b are real numbers characterising the action. It also acts on the complex 
numbers: 

s:R + xC^C, s(n,z) = s^z) = fi c z, (13) 

where c is yet another real number. Now, since the group acts on M 4 (via a) and 
C (via s) it also acts on I 4 x C (via a x s) and hence on graphs of functions 
^ : M 4 — > C. The transform of a function (via S) is then defined to be the 
function uniqueley corresponding to the transformed graph: 

graph(S /J *) : = cr M x s M (graph(*)) 

(t,z)eR 4 (14) 

(t,x) 



(15) 



where in the last step we just used that each is a bijection. Hence 

S : M+ x FunctQR 4 , C) Funct(M 4 , C) , 

%$) = S,($) = Stl o*o ff( ; 1 . 

or, in simpler terms, 

S^(t, x) = i^ c ^(fi- a t, pT b x) . (16) 

As already mentioned, we wish to consider normalised wave functions. The 
scaling ( 15 ) preserves normalisation iff 2b + 3c = 0, i. e. iff c = —36/2. Solu- 
tions to the SN-equation for all positive mass parameters m then correspond to 
special paths on the unit sphere in Hilbert space. We ask under what conditions 
the action ( [T3] ) transforms solutions ^ to the SN-equation for mass parameter m 
to solutions S^ty for some mass parameter m M . We have 



dtSpV = ^ a S^ (17a) 
m~ 1 AS r M * = ^ 2b {m/m^) m^S^A® (17b) 
mlU(S^)S^ = (mjmf m 2 S^U (17c) 
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where U(&) is the potential function 



U(*)(t,£) = I ^4fd 3 x. (18) 



\x - y\ 

Denoting the SN-equation symbolically by SN(m, \I/) = 0, we have 

SN(m,^) = => SN(m^S^) = (19) 



iff each term scales with the same factor. According to ( 17 1 this is the case iff 
H~ a = ii- 2b (m/m M ) = iT h {mjm) 2 (20) 

i. e. iff 

^ = /x -*>/3 = Ar a/5. (21) 
m 

First of all this means that is a power of /i times m. This power may without 
loss of generality be taken to unity, for if m^/m were equal to {i d we could take 
fjf := fJL as new and equally good scaling parameter and everythig that follows 



were true for // instead of /i. So from m^/m = fj, we see from (21 ) that b = — 3, 
a = —5, and hence c = 9/2. To sum up, we have shown that 

S^(t,x) = ^'H{iih, fi 3 x) (22) 

satisfies the SN-equation for mass = [i-mif^ satisfies the SN-equation for 
mass m0 

The stationary version of ([[]) is obtained by replacing ihdt with the energy E. 
It has a unique (up to translations, of course) spherically symmetric stable ground 
state. The ground state minimises the energy functional which is given by the negative 
of the second and third term in Q. Existence and uniqueness of the ground state 
were first shown by Lieb |6], albeit motivated in completely different physical context 
(one-component plasmas) in 1976 by Choquard. Hence Lieb called it the Choquard 
equation, which is mathematically equivalent to the stationary SN-equation. 

Numerical studies led to the following value for the ground-state energ>|^] 

G^TfC* ( Tfi \ 

E Q = -0.163— j— = -0.163 -mc 2 • , (23) 

n z \mp J 

where 

[he- 
mp ■= J = 1.221 x 10 19 GeV/c 2 = 21.76 /zg (24) 
V Cx 

is the Planck mass. A rough intuition for the width, a, of the mass distribution m-l^l 2 
in the ground state is obtained by schematically writing the total energy as sum of 



Formulae (2.3) of |9| and (2.4) of [2] seem to be incorrect. 

See |8| for an early attempt and |5| for a more recent and far more detailed account, also showing 
linear stability of the lowest and linear instability of all higher spherically symmetric stationary states 
found earlier in [7 |. The 20 lowest energies are tabulated in Table 1 of | 5 1 and range from —0.163 to 
-0.000221 times G 2 m 5 /h 2 . 
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the kinetic and potential term and replacing the Laplacian in the first by a 2 and the 
double integral in the third term of (Wj) by a -1 : 



(25) 

2ma 2 2a 



«o«^ = 2^p-( — ] • (26) 



Minimising in a then gives 

Gm 3 v m 

where 

£ p .-J = 1.616 x 1(T 35 m (27) 



is the Planck length. The rough energy ( 25 ) at the rough minimum ( 26 1 is 

E °~-8'^- (28) 



which compares reasonably with (23 ). 

Clearly, we would not trust the stationary solution in this Newtonian picture if the 
support of the mass distribution m| , I'| 2 had a spatial extent not much larger than m's 
Schwarzschild radius Gm/c 2 . Now, the condition ao 3> Gm/c 2 is equivalent to 

m 4 < 2 • m% , (29) 

which means that we should not let m come too close to the Planck mass (though note 
the 4th powers). 

Finally we note the dimensionless form of ([TJ that one obtains by introducing a 
spatial length scale I. Using dimensionless spatial coordinates x' and a dimensionless 
time t' according to 



VI, (30) 

and rescaling ^ h- > Vl/' := i^l 2 ■ ^ so as to keep ty' normalised with respect to d 3 x', 
we get 

ifyVV,*) = f-A' - K J K^S dV) *'(t',x') (31) 
with dimensionless coupling constant 

Note that the combination (length) x (mass) 3 which determines the coupling had to 
be expected from scale invariance ([22b. 
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3 Analytical Estimations of the Collapse Mass 



We consider now the time-dependent SN-equation for initial values given by a spheri- 
cally symmetric Gaussian wave packet of width a: 



*( r , t = 0) = (Tra 2 )" 3 / 4 exp U^J . (33) 

We have then two free parameters, a and m, and we ask for the regions in this two- 
parameters space in which significant inhibitions of the usual free dispersion occur. 
For simplicity we will refer to such a behaviour as a "collapse", even though it is 
clear that neither the integal over |^| 2 nor the volume of the support of \ip\ 2 (which 
is infinite) change. We start by presenting four different analytical arguments, all of 
which suggest that a collapse of the wave packet sets in whenever the dimensionless 



coupling K in ( 32 1 approaches unity. 



3.1 Equilibrium of Forces for the Free Solution 

Let us first give a rough estimate for the critical mass by a naive calculation, in which 
we demand equilibrium between the acceleration of the peak probability density and 
the acceleration due to gravitation for a point mass. 

The free solution (without gravitational potential) for the problem at hand is 

WM) = (™r 3/1 (l + S)" 3 ' 2 exp (- 2a2(l "^ ) ) ■ (34. 
The radial probability density, p(r, t) = An r 2 l^ftee^, t)\ 2 , has a global maximum at 



\ 1 • ''Tr (35) 
Taking the second time derivative we obtain the peak's acceleration 

h 2 / h 2 t 2 y 3/2 h 2 

m z or \ m z ar J m z r£ 

This we now compare to Newton's law for the gravitational acceleration 

Gm 



(37) 



at time t = 0. Using r p (t = 0) = a, we get 



» 2 ^ 1/3 



m= [a-a) ■ <38) 

For a = 0.5 /jm this yields a threshold mass for collapse of about 7 x 10~ 18 kg or 
4 x 10 9 u. 
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3.2 Lower Mass Bound for Negative Energy Solutions 

It was noted by Harrison et al. [5 ] that only initial data of negative energy can show 
inhibitions of dispersion. This can be used to provide a lower mass bound for collaps- 
ing behaviour in a mathematically rigorous way. The result strongly contradicts the 
collapse mass found in Sill- 

Define the kinetic energy T, the potential energy V, and the total energy £ respec- 
tively by 



£=t + \v 
h 2 



2m 



j d 3 x|W(x)| 2 + y J d 3 x<S>\^(a 



, d 3 x 2 / d 3 x / d 3 ?/ 1 V-! -V ■ 

2m ./ 2 7 J x-y 



The second moment, defined by 



(39a) 



Q 



\x\ 2 |^| 2 d 3 x, 



then has the following second order time derivative: 



Q 



2h 2 



m c 



|W| 2 d 3 x + / $|tf| 2 d 3 x 



1 



m 



(4T + V) 



m 



(40) 



(4£ - V). (41) 



But (and therefore also V) is everywhere negative due to the maximum principle 
which means that for positive energy £, Q is a convex function of time. Hence, in terms 
of the second moment it is a necessary (but not sufficient) condition for a collapsing 
wave packet, that its initial data has negative energy. 

We can now easily calculate the energy of the spherically symmetric, real Gaussian 
wave packet (|33]>. It is 



£ 



r h 2 f°° 

/ dr(d r ^(r)) 2 -8ir 2 Gm 2 

n Jo 

f-oo r rr 12 roo 

dr dr' — ^ 2 (r)^ 2 (r') + / dr'r'f 2 (r)f 2 (r') 

J UO r Jr 



2h 2 



0rma 7 J 



drr 2 e -r 2 /a 2 



8Gm 2 



dr 



-r 2 /a 2 



/ j / 12 -r' 2 /a 2 i / j / / 

dr r e ' + r dr r i 

JO Jr 



-r' 2 /a 2 



h 2 



2Gm 2 



2ma 4 



ir or 



Now £ < is equivalent to 



m > 



sinh _1 (l). 



n 2 



(42) 



AG a sinh-Vl] 



1/3 



ft 2 \ 1/3 



2Ga 



(43) 



By inserting a = 0.5 again this yields the lower mass bound for negative energy 
solutions m > 5.5 x 1CT 18 kg « 3.3 x 10 9 u. 
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3.3 Mass of the Stationary Ground State 

In our third example we consider the ground state of the spherically symmetric SN- 
equation found by Moroz et al. in Q. In order to do this we need to rescale the 
SN-equation in a suitable fashion in order to to compare our form of the equation, 
which explicitly contains the mass parameter, with the form solved in Q, in which m 
does not appear. In order to understand how to transform solutions of one form into 
solutions of the other, we use scale invarianvce. 

In a first step we substitute U = j£(E — m<&) into the stationary SN-equation 



h 2 



2m 



E^(x) 
A$(z) 

It then takes the following form: 

= -M>(x) , 
2m 2 



A + m&(x) 



47rGm|^(f)| 2 



AU(x) 



h 2 



8vrGm 3 

h 2 



This can be further simplyfied if written in terms of the rescaled function Q 
where a := V 8irGm 3 /h. We get 



AO(f) 
AU(x) 



-U(x) n(x) , 
-\V(x)\ 2 . 



(44a) 
(44b) 

(45a) 
(45b) 

= a®, 

(46a) 
(46b) 



The norm of Q is obviously = a \\^f\\. 

Now we rescale the mass by a factor of ji, as discussed above. Then, according to 



equation p2), S^(x) = n 9/2 ^f(fi 3 x) and a 



a. 



3 / 2 <x Therefore 



(47) 



n 6 n(n 3 x), 
M 3/2 ||n||. 

This means that the rescaling changes the norm of Q and the following statement holds: 



Sfj,Q(x) 

IIV^II 



/i 6 a^(/x 3 f) 
/i 3/2 a H^ll 



Proposition 1. Suppose that (0, U) satisfy ( |46| ). For a given mass value m let ji be 
the rescaling parameter for which the norm of S^Q is 



Then the wave function defined by 

*fx) : 



S^ix) 



(48) 



(49) 



is a normalised solution to the stationary SN-equation ( 44 ) with mass parameter m. 



Putting equations (|47|)-(|49|> together we obtain a formula that gives us a normalised 
wave function satisfying the stationary SN-equation ( 45 1 with mass parameter m for 
any unnormalised solution Q of ([46]): 



(8ttG) 3 / 2 m 9 / 2 (SnGm 



ft 3 imn 4 



n 



h 2 \\n\\ 2 



(50) 
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O(r) 
1.1 r 
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Figure 1: The ground state solution (|5l} with Q = 1.08864 found by Moroz et al. [7] 
plotted up to forth order. It seems reasonable from this plot to cut off the solution for 
values of r > 2.2. 



In Q Moroz etal. found spherically symmetric, real solution^] 

1 



n(r) = n i 



_ r 2 , ^0 + 1 r 4 

6 120 



(51) 



to (45), where f^o = 0(0) and £Iq rs 1.08864 for the ground state. These polynomial 



functions are in general not normalisable but for the right choice of Uq they tend to 
zero asymptotically as r — > oo. Nevertheless, for the approximate solution we must 
introduce a cutoff A such that Sl(r) = for r > A. The ground state solution is plotted 
in figure[T]up to order r 4 , according to which A 2.2 seems to be a reasonable choice. 
Therefore 



\n\\ 2 » [ 47rr 2 drft 2 (r) ~ 26. 
Jo 



If we define as the width of Q the cutoff A, the width of the wave packet ^ is 



h 2 \\n\ 



A 



57 h 2 



or written in terms of m 



87rGm 3 8vrGm 3 ' 
57h 2 \ 1/S 



m 



8nGa 



With the usual value of a = 0.5 we get a mass of 



m 



9.1 x 10~ 18 kg ps 5.5 x 10 y u. 



(52) 



(53) 



(54) 



(55) 



For the spherically symmetric SN-equation, the solution can w.l.o.g. be assumed real. This follows 
from J45al if restricted to the spherically symmetric case. 
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3.4 Short-Time Behaviour 



Finally, to obtain the behaviour of the spherically symmetric SN-equation with an 
initial Gaussian wave packet for small times, we consider the lowest order terms of the 
series expansion around t = 0. Transforming to (f>(r,t) := r^(r,t) and considering 
the spherically symmetric gravitational potential 



$(r,t) 



-AirGm 



- [ W,t)\ 2 dr' 
r Jo 



+ 



(r',i)|< 



dr' 



the SN-equation together with its complex conjugate reads 

ih 



d. 



ih 2 im 



i m 



(56) 



(57) 



where, as before, the star denotes complex conjugation. We get the first order time 
derivative 



ih 



(58) 



of the absolute value squared which vanishes at t = because <fi(r, 0) is real valued. 
For the second order time derivatives we obtain 



''' ;, 

2m 



— (d t ®) <(> ~ -y $ (dt<f>) 



h 2 



1 



m 2 2 



2 /* 



d 2 \6\ 2 



h 2 



1 



• m2 + g (^<*>)<f + (WW) + * (3 



2 .* 



(^)f + 2(W) + 

"^2 0* " 2 (5 r ^) W) + W)] 

+(d r 2 <I>) |</>| 2 + (W ( W 0* + (9 r ^*)) . 



(59a) 



(59b) 



(59c) 



We introduce some shorthand notations 



a 
71 



(vra 2 )^/ 2 
4vr- 1 / 2 Ga/3 
—a mT 1 $ 



and consider the value of d 2 \(j)\ 2 



P:=a- 1 p:=/3r 
7 2 :=2Ga/3 73 := h 2 a (3 2 
^o(r) :=0(r,O) = ^re"^ 2 / 2 , 

at t = which, as $ ~ m, takes the form 
1 



9t\4>\ Z \ t=0 = mA(p) + -^B(p), 



(60a) 
(60b) 
(60c) 



(61) 
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(a) A(p) 



(b) B(p) 



Figure 2: Components A(p) and B(p) of the order i 2 contribution to \<f>\ 2 which is pro- 
portional to m A(p) + -K B(p). While A is attractive, B shows a repulsive behaviour. 



where 



A(p) = ((<9 2 <D) 0g + 2 (a r c&) (d r cp ) O ) (62) 
#(p) = -y (<Ao (5 r Vo) - (<9 2 0o) 2 ) . (63) 



Evaluating the derivatives of $ we obtain 





$ = 


i r , 2 - 

- / p e 
P Jo 


/•oo 

J p 






erf(p) 








4 p 






= a p $ = 


2p 


spK erf(p) 
4 p 2 


^* 


= d 2 $ = 


\pK erf(p) 
2 p 3 





where erf (p) is the error function 



o 



Inserting (64b I and (64c I into (62i, the final result can now be stated as 



(64a) 
(64b) 
(64c) 



erf(p) := ^= /' e'^dx. (65) 



H 2 = </>o + ^ (mAW + iB( P )j+0(t 3 ), (66) 

where 

A(p) = 2 7i p 2 e" 2 " 2 - 72 p e"^ erf (p) (67) 

B(p) = 73P 2 (2p 2 -3)e-" 2 . (68) 

A(p) and B(p) are plotted in figure]!] While A shows an attractive behaviour (i. e. the 
wave function grows for small values of p), B shows a repulsive behaviour. Thus, the 
attractive behaviour gets dominant for large masses, as expected. 

In figure [5] the second order derivative of \<p\ 2 is plotted for several values of the 
mass m. One can see that there is a repulsive behaviour for m < 4 x 10 9 u which 
turns with increasing mass and is completely attractive for masses m > 9 x 10 9 u. We 
therefore expect the collapse of the wave function to happen at about 5-8 x 10 9 u. 
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(a) m = 10 b u 



(b) m = 4 x 10 9 u 



1*1- a, 3 i0i ; 




(g) m = 9 x 10 9 u (h) m = 10 12 u 

Figure 3: Plot of <9 t 2 |</>| 2 at t = for several masses. While the graph changes only in its 
quantitative scale and does not change qualtitatively in its shape for masses smaller than 
m = 4 x 10 9 u and larger than m — 9 x 10 9 u, respectively, one can see a qualitative 
change from repulsive to attractive behaviour in between those mass values. 
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Figure 4: Collapsing wave packet for m — 7 x 10 9 u. Plotted is the radial probability 
density p — 4tt r 2 \ip\ 2 against r at three different times. 



4 Numerical Investigation of the time dependent Schrodin- 
ger-Newton Equation 

Our numerical methods do not differ in principle from those used in J5J. A Crank- 
Nicolson scheme is applied to solve the differential equation ([T]) in a spherically sym- 
metric context. 

Having defined both a spatial and temporal grid size Ar and At and using the 
index notation \I>" = ty(jAr, nAt) we can write equation ([T]) in a discretised way 
using Cayley's form of the SN-equation 

i At 

exp 



Linearising this equation we can write it as 

Q 



2h 



H) *?. 



2 V 2ft 



We thus have to solve the linear system 

Qx n = * n 

to obtain 

The radial component of the Laplacian in spherical coordinates, 

r "\3^ ifr-0, 



(69) 

(70a) 
(70b) 

(71) 

(72) 

(73) 
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Figure 5: The peak of the radial probability density plotted against time for several 
masses. The curves are marked with the mass in units of 10 9 u. 



takes the discretised form 



(Arp 



while the discretised form of the potential 

1 



' \) i ^ • 7 \].) if J >0 



if J = 0, 



is 



<j> n 

3 



\n(Ui, | = J \^(r',t)\ 2 r' 2 dr' + 



|*(r',t)|V dr 



-ATiGm{Ar) 2 v] 



iV-l 



Q then becomes a tridiagonal matrix 





(bo 


co 





••• \ 
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Figure 6: The radius rgo within which 90 % of the probability lie plotted against time for 
a mass of m — 10 10 u. Note that the minimum is not at zero but at about rgo ~ 0.4 /im. 
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We use the tridiagonal matrix algorithm to solve the linear system (|7T 
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4.1 Numerical Results 



We start with a spherically symmetric Gaussian wave packet (33 1 whose width is al 



ways set to a value of a = 0.5 /xm. Note that due to the scaling law ( [22] ) the width can 
be fixed without loss of generality. 

The results can be summarised as follows: 

• For masses m < 6 x 10 9 u the wave packet spreads, but slower than the free 
solution for zero potential spreads. The larger the mass, the slower the spreading 
becomes compared to the free solution. 

For masses below 10 9 u the difference to the free solution is hardly identifiable. 

• For masses m > 7 x 10 9 u we observe a collapsing wave packet. The shrinking 
behaviour was verified for masses up to 3 x 10 11 u. 

• We have no results for larger masses, yet, as for reliable results one needs to in- 
crease both the temporal and spatial resolution with increasing mass, and there- 
fore the runtime of the calculation also increases. 

In figure |4] snapshots of a typical collapse are presented. The radial probability density 
p = 4-7T r 2 |^| 2 is plotted against r for a mass of m = 7 x 10 9 u. 

In figure [5] we plot the time evolution for the global maximum of the probability 
density. The collapsing behaviour for masses not less than 7 x 10 9 u is obvious. Also 
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Figure 7: Time it takes until the gravitationally interacting solution differs from the so- 
lution of the free Schrodinger equation in its FWHM by a percentage of 1 %, 10%, and 
50 % respectively. 
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note the oscillations for larger time values which are due to the packet penetrating 
itself. There also seems to be some equilibrium distance to which the wave packet 
collapses, which might correspond to a stationary solution. The chaotic behaviour for 
the curve corresponding to m = 5 x 10 10 u is probably due to numeric effects. 
We also plot the radius rgo within which 90 % of the probability are located, i. e. 

90%= / 47rr 2 |^| 2 dr, (79) 
Jo 

against time. One can see from figure |6j where rgo(i) is plotted for a mass of 10 10 u, 
that it is reasonable to define a collapse time by the minimum of this plot, which then 
is about 30 000 seconds. 



4.2 Coherence Time 

While interferometry experiments might still be feasible at masses of order of the 
collapse mass of 10 10 u that we found, another experimental constriction occurs. As 
can be seen from figures[5]and|6]the time scale beneath which the collapse takes place 
is of the order of several hours. For experimental tests one would like to take this time 
scale down to a reasonable value. 

Hence, we plot in figure [7] the time it takes the solution of the full SN-equation to 
differ from the solution of the free Schrodinger equation in its full width at half max- 
imum by a certain percentage against the mass. Figure |7](a) shows that the coherence 
time needed decreases with increasing mass. According to figure |7](b) it takes, for 
example, about 3 minutes for the solution for m = 3 x 10 11 u to differ from the free 
solution by one per cent. 

Consequently on one hand, if keeping the width of the wave packet constant, for 
an experimental test of the SN-equation by means of e. g. molecular inteferometry 
one needs to maximise not only the mass but also coherence time and the sensitivity 
regarding the detection of deviations from the free Schrodinger equation. 



But note that on the other hand, we also could make use of the scaling law ( |22| ) 
to decrease the collapse time. If we simultaneously increase the mass by factor of n 
and decrease the width by a factor of the collapse time decreases by a factor of 
/i -5 . Thus, for a mass of 10 11 u, for example, observed with a grating period of 0.5 nm 
one should in principle be able to observe a collapsing wave packet causing a loss of 
interference with a coherence time of approximately 300 ms. 
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